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Mafia (also called Werewolf) is a party game. The participants are divided into two competing 
groups: citizens and a mafia. The objective is to eliminate the opponent group. The game consists 
of two consecutive phases (day and night) and a certain set of actions (e.g. lynching during day). 
The mafia members have additional powers (knowing each other, killing during night) whereas the 
citizens are more numerous. 

We propose a simple mathematical model of the game, which is essentially a pure death process 
' with discrete time. We find the closed- form solutions for the mafia winning-chance w(n,m) as well 

as for the evolution of the game. Moreover, we investigate the discrete properties of results, as well 
as their continuous-time approximations. 

It turns out that a relatively small number of the mafia members, i.e. proportional to the square 
root of the total number of players, gives equal winning-chance for both groups. Furthermore, the 
C^) • game strongly depends on the parity of the total number of players. 

Oh 

^3 ■ i. introduction 

Mafia (also called Werewolf) is a popular party game The participants collect in a circle and a game 
coordinator assigns each player to one of two groups: a mafia or citizens. Citizens know only their own 
, identity, whereas mafia members know identity of their fellows. The objective is to eliminate the opponent 

£> ' group. The game consists of two alternating phases (day and night). During the day, all players can discuss 

and vote who they want to lynch. During the night the mafia kills a citizen of their choice. 

cn . 

There are many variants of the mafia game. The most common modification of the gameplay is the addition 
of characters with special abilities. Usually there are two special citizens: Detective (who checks every night 
, if a chosen person is in the mafia) and Nurse (who may protect a victim from being killed, if she chooses 

f*^ ■ correctly). 

In this paper we analyze the simplest version of the Mafia game. That is, with only mafia and citizens, and 
without any special players or additional rules. Even if it may not be the most popular variant, it is the best 
one for mathematical modelling. Despite the game having a complex psychological component it is possible 
to create a stochastic model of the Mafia game. We consider a game beginning with n players, out of which 
t , ■ there are m mafia members. The main questions we address are: 



o 



• What is the probability that the mafia wins w(n, m] 



? 



• How does the dynamics of the game look like? That is, what is the chance that after a given time 
there is exactly a certain number of mafia members? 

Besides the direct answers (i.e. closed form expressions), we study approximations, qualitative behavior and 
some special cases. 

There are only a few previous papers on the Mafia game. Works [2| and [3| basically propose the same 
model as the one presented in this paper. They calculate a simple asymptotic formula for the mafia-winning 
chance, w{n,m) oc m/^/n. 



'Current address: ICFO— Institut de Ciencies Fotoniques, Av. Carl Friedrich Gauss, num. 3, 08860 Castelldefels (Barcelona), 
Spain; Electronic address: pmigdal@gmail.com, URL: http://migdal.wikidot.com/en 
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Other research projects concentrate on the psychological aspect of the Mafia game, in particular — deceiving. 
Topics analyzed include: patterns of a deceiver's interruptions and voice parameters 1 , movement of face 
and hands [5| and usage of language [6]. Also, the Mafia game is applied as a test for some videoconferencing 
setups Q. 

The research is motivated by two goals, besides the sheer fun of calculation and writing. First, a good theo- 
retical model may be useful for the investigation of the psychological aspect of the Mafia game. Experimental 
deviations from the idealized behavior may give a valuable insight into psychology of strategy choosing, ma- 
nipulation, deceiving, following others and hiding identity. Second, when a small, but well-informed and 
powerful group fights against the majority, the nature of the process may be similar to that of the Mafia 
game. Such real-life phenomena include actions of secret societies as in high-stake corruption, terrorism and 
illegal oppositions. 

The paper is organized as follows. In Section [TT] we introduce the rules of the Mafia game and propose a 
simple mathematical model. Section IIIII contains a special case of play with only a single mafia member. 
The results are not only simple and didactic, but also useful for a more general case. Section IIVI gives 
simple results on the qualitative behavior of the mafia- winning chance. Section [V] presents the dynamical 
aspect of the Mafia game and gives closed-form expressions for dynamics of the game. We consider both the 
discrete-time model and its continuous time approximation. In Section I VII we provide the exact result for 
the mafia- winning chance, together with its asymptotic approximation. Section [VIII concludes the work and 
gives insight into possible extensions and applications. Appendices |A] and [B] contain derivations of formulas 
found in Sec. [V] 



II. MODEL 

Before proceeding to mathematics, we need to write down our arbitrarily chosen rules of the Mafia game. 
Some assumptions are made to simplify the model, others purely for convenience. Note that this Section is 
only a hand-waving transition from the real-life game to its mathematical model. 

• The game needs n players and one more person to coordinate it. 

• At the beginning players are randomly divided into (n — to) citizens and to mafia members. 

• Mafia members know the identity of each other, citizens — only their own. 

• There are two alternating phases, day and night, which together comprise a turn. 

• During day, there are two consecutive subphases: 

— Debate. Everyone still alive can say anything related to accusing or defending. 

— Voting. Everyone has exactly one vote, who he wants to lynch. The player who gets the highest 
number of votes is eliminated (in case of a tie, a random 'winner' is eliminated). 

• During the night: 

— Mafia jointly decides who they want to kill (eliminate). 

• The game continues until there is only one group (either the citizens or the mafia) left. That group 
wins. 

Stated briefly, during the day one player is eliminated (either a citizen or a mafia member) whereas during 
the night the mafia kills one citizen. Consequently, during a single turn the possible transitions are: (n, to) — > 
(n — 2, to) and (n, m) — > (n — 2, m — I). Their probabilities, in principle, depend on n, to and the course of 
the play (i.e. previous discussions and votings). To start with, we need to get rid of the psychological aspect, 
restricting ourselves only to strategic and probabilistic parts of the Mafia game. It is a coarse approximation, 
as psychology plays an essential role in the Mafia game. Nevertheless, even such a bare model has interesting 
properties. Moreover, we want to neglect the discussion phase as it is extremely difficult to formalize in a 
meaningful and useful way. 
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Now, let's make a hand-waving argument why the transition probability depends only on the current state 
(n,m). Consider the scenario, in which players are lynched at random, that is 

77 — 777 777 

P\(n- 2,m)\(n,m)] = and P \(n - 2, m - l)\(n, m)l = — . (1) 

77. 77, 

Let's denote by w(n,m) the mafia-winning chance, induced by the above transition probabilities. 

Suppose that there is a better strategy (in the common sense of the word) for citizens, giving it w'{n, m) < 
w(n,m). Then the mafia can force a random lynch, just by pretending they are citizens. It is possible as: 

• Each mafia member knows more than any citizen. 

• In such case there is no way to detect a mafia member. 

Suppose that there is a better strategy for the mafia, w"(n,m) > w{n,m). If the citizens are in majority, 
they can force a random lynching by the following procedure: 

• One citizen says "let's every of us give a number ki from 1 to 77, we will kill the ^2 ki mod 77-th". 

• The number ^2 ki mod n is random, so is the chosen player. 

• As there are more citizens than mafia members, mafia votes do not play any role. 

If the citizens are in the minority (i.e. m > n — m) no matter what the moves are, the mafia will win. For 
777 = 77 — 77i let's assume a random player is lynched. 

Consequently, random lynch is a well-justified strategy. A more detailed argument is shown in [2]. As a 
direct consequence of ([1]), the mafia winning-chance can be expressed as a recurrence equation 

if m = 0, 

if m > n — m, (2) 
2, m) + —w(n — 2, 777 — 1) in all other cases. 

The above equation is the core of our paper. We investigate it for one mafia member (Section IIII[) . its 
qualitative behavior f Section IIV|) and closed- form solution along with its asymptotic formula (Section IVI[) . 

It is important to point out that we may set different boundary conditions, i.e. in which if there is the same 
number of citizens and mafia members, mafia wins. That is, in @ we may write (1 if 777 > 77 — m) instead 
of (1 if to > 77 — to). In that case all results can be reproduced, perhaps in a slightly changed form (i.e. 
different constants and evenoodd). 





Ui(t7,, 777) = < 1 



-w{n — 




Figure 1: A numerical plot of the probability that the mafia wins w(n,m). 
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III. GAME WITH ONE MAFIA MEMBER 



To start with, we want to consider a very simple case — the game with only a single mafia member. That 
is, in this Section we analyze w(n, 1). Besides the simplicity, there are two more motivations: this case gives 
insight into some general properties and the results will be useful in the further part of the paper. 

In a game with one mafia member, their winning probability is the chance that during every lynch a citizen 
is killed. Bearing in mind that we start the game with the day phase, we get 

n-1 rc -3 1 + (nmod 2) _ (n- 1)!! 

______ _ . (3 ) 

The above formula has an explicit dependence on the parity of n. Even though it is obvious that addition of 2 
citizens reduces mafia-winning chance w(n + 2, 1) < w(n, 1), it is not the case for addition of a single citizen. 
In fact, addition of an odd player (holding m — 1) increases the mafia- winning chance. To understand this 
let's use an example. In a play with one citizen and one mafia member, there is tie in voting, so lynching 
relies on tossing a coin and w(2, 1) = 1/2. In a game with two citizen, the mafia also wins after lynching a 
citizen (the second one is to be killed during the night). This time it is more difficult to hunt mafia, and 
w(3, 1) = 2/3. One may think that we encountered an issue of boundary conditions (i.e. killing a random 
player in a tie) or that we considered too small n. Neither is the case. A careful reader may easily check 
that changing (1 if m > n — m) to (1 if m > n — m) in ^ we get an analogous phenomena. 

We may simply check that adding an odd player always increases the mafia winning-chance, 

, , x , , x 2fc 2k - 2 4 2 

w(2k + 1, 1) - w(2k, 1) = ... + 4 

v ' ; v ' ' 2k + 1 2k- 1 5 3 1 ' 

2k -I 2k -3 3 1 



2k 2k -2 4 2 



v/(2fc + l)(2k- 1) y/(2k-l){2k-3) V5~3 V3~l 
> 2k + 1 2k- 1 '"' _ 5 3 

2k - 1 2k - 3 3 1 



^{2k + l)(2fc- 1) ^/(2fc — 1)(2A; — 3) 



= 0, 



that is, w(2k + 1, 1) > w(2k, 1). 

Let's consider w[n, 1) averaged over neighbouring numbers. We use the geometric mean, as it simplifies the 
result 



i—, — r / (n - 2)!! (n - 1)!! 1 



(5) 



The averaging may be considered an approximation in which we get a monotonic function (without the 
sawtooth pattern). To get some information on the dependence of w(n, 1) on the parity of n we consider the 
following ratio 

w(2k + l,l) = (2k)\\ (2fc-l)!! = / (2fc)H \ 2 1 
_(2fc,l) (2fc+l)!! 7 (2k)\\ \{2k-l)\\) 2k + V U 

The above expression has a limit when k goes to infinity, which can be found with the aid of the Wallis 
formula: 

7T / (2fc)H \ 2 1 

2 ~ fc^o ^(2fc- 1)!! / 2fc + 1' U 



We may develop approximate formulas for the single mafia member winning- chance, which take into account 
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parity of the number of players 



io(2fc,l) 



Or, bearing in mind (n + 1) /n — > 1, 



^) 4 v^WiHi) (8) 



2 1 



7T 1 



2 V2fc + r 



/7r Wnmod2)-l/2 1 



As we see, the recurrence equation produces results that might be counter-intuitive. Imagine there are four 
of us and we are playing with a single mafia member. His winning-chance is iu(4, 1) = | = 0.375. Then 
we invite five more friends to play as additional citizens. Consequently, the mafia winning-chance rise to 
w(9, 1) = I?! « 0.406 (contrary to a naive expectation). 



w(n,l) 
1.0f-.\ 
0.8 [ K 
0.6 
0.4 
0.2 

5 10 15 20 n 




Figure 2: Plot of the mafia winning-chance for a game with a single mafia member, that is w(n, 1). Dots illustrate 
the exact result ([3]), whereas the lines represent the approximations ((9} — solid and dashed line for even and odd 
number of players, respectively. 



IV. QUALITATIVE PROPERTIES OF w(n,m) 

It is worth investigating the qualitative behavior of w(n, m). Some questions, which naturally arise, are: 

• Does the addition of two citizens increase their winning-chance, i.e. is w(n + 2, m) < w(n, m)? 

• Does the addition of one citizen and one mafia member increase the mafia winning-chance, i.e. is 
w(n + 1, m + 1) > w(n, to)? 

• Does the change of one player from a mafia member to a citizen decrease the mafia winning-chance, 
i.e. is w(n, to + 1) > w(n, to)? 

In fact, the above questions are equivalent, as we can see using using the recurrence relation @ 

. n — to . to 
w(n, to) = w(n — 2, m) H w(n — 2, to — 1) 



G 



for n — to > to > 0. After some cosmetic arithmetical operations we get 

(71 — 771) [w (n — 2, m) — w{n, m)] = m [w(n, m) — w(n, m — 1)] , 

ri [tx>(71, 771) — 7x1(71 — 2, m — 1)] = (71 — m) [w{n — 2, m) — w(n — 2, 771 — 1)] . 
The above relations can be directly translated into 

tx(?i — 2, 771) > w(n, 771) 7x1(71, 771) > tt;(7i — 2, 771 — 1) <=> w(n — 2, 771) > io(n — 2, m 



(10) 



(11) 



In addition to the equivalence of the 'obvious' inequalities, we are interested whether they are fulfilled. Let's 
utilize the third part of (fTTj) in the following form: 

If a mafia member is changed into a citizen, the mafia winning-chance does not increase, that is ix>(ti + f , tti + 
f ) > w(n + 1, 771). If additionally n — m > m, the inequality is strong, 7x(7i + 1, 771 + 1) > 7x>(ti + 1, 771). 

Let's introduce an agent to a game with 71 players among which there are m mafia members. The agent is 
an additional player who simulates either a citizen or a mafia member. We are going to show that the agent 
does not have to reveal its identity until the very late part of the game. 

As citizens are indistinguishable, let's assume that mafia kills them during night in a fixed sequence. Let 
the agent have the last position. During the lynching there are three possibilities 

• A citizen dies. 

• A mafia member dies. 

• The agent dies. 

In the case of the agent's death, his affiliation has no meaning. In the two remaining situations, the game 
continues with the agent. The agent's affiliation plays a decisive role in two situations: 

• When all mafia members are dead. 

• When 0-1 citizens and 1-2 mafia members remain. 

In both cases when the agent is a mafia member, the mafia winning-chance increase: 



Game state (77, 771) 


7x> if the agent is of mafia 


w if the agent is of citizens 


(71,0) 





?x(7i, 1) 


(1,1) 


1 

2 


1 


(2,1) 


1 


1 


(3,2) 


1 

4 


1 



The total mafia winning-chance is equal to the mean winning-chance, averaged over disjoint games. The 
weight is the probability of achieving a such game. The change of the agent's affiliation does not affect 
the weights, but in the ending always increases the mafia winning-chance, 7x>(?i + 1,777, + 1) > w(n + l,m). 
Additionally, for 71 — 771 > 771 > 1 there is a game in which all mafia members are killed and the agent's 
affiliation plays a role, hence w(n + 1, m + 1) > tx(?i + 1, 771). 

Let's show one more property — when 71 — 771 > 771 adding an odd player always increases the mafia winning- 
chance, not only in games with one mafia member (|4]). We will use mathematical induction, using the 
boundary conditions ([2]) and the already shown properties of 7x;(n, 1) as the basis. The inductive assumption 
is that the property holds for the total number of players n < 2k + 1, and the inductive step is the following: 



;(2(* + l) + I,™) = 



> 



> 



2fc + 3 


— 771 


2k + 


3 


2k + 3 


— 771 


2k + 


3 


2k + 2 


— 771 


2k + 


2 


w(2(k- 


hi), 



?x;(2fc+ l,77i) + 



2k + 3 



7X>(2fc+ 1,771-1) 



(12) 



7X>(2fc, 77l) 

w{2k, m) 
771), 



2fc + 3 

771 

2fc + 2 



w(2k, 771—1) 
w{2k, 771—1) 
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where we subsequently used the inductive assumption and made use of the fact w(2k, to) > w(2k, m — 1). 



V. EVOLUTION 



In the two previous Sections we considered the mafia winning-chance with respect to the initial number of 
players (N) and mafia members (M). It may be interesting to explore dynamics of the play, that is, analyze 
the probability p m (t) that after the t-th turn there will be exactly to mafia members. Even though this goes 
beyond the ([2]) , it is a direct consequence of the random lynching model (fTJ) . Throughout this section we use 
capital letters N and M to denote the initial conditions. 

Note that every turn, or day and night, two players are killed, so the total number of player decreases 
n(t) = N - 2t. 



A. Discrete time 



The evolution of each p m (t) is governed by a set of recurrent equations 

(N-2t)-m , s m + 1 , , .„„, 

Pm(t + I) = N _' 2t Pm(t) + Jj—^ Pm+1 (t), (13) 

where to € N and with the initial condition p m (0) = S m M- 

Note that the above is a stochastic process with M + 1 distinct states and (to) — > (to — I). In general, such 
a Markov process is called a pure death process (a subclass of birth and death processes) . In our case the 
transition probabilities change with time, which is a slight complication. 

It is essential to point that the equations in (jI3[) are correct only when m < n{t), that is, when there is at 
least one citizen alive and the game can be still played. Otherwise they have no meaning. But why does 
an erroneous equation for m > n(t) not spoil states with m < n(t)7 For a given time t let p m (t) be the 
erroneous probability with the smallest index, i.e. m = N — 2t + 1. In the next turn it affects p m _i(t + I). 
But, now the smallest wrong index is N — 2(t + l) + l = m— 2 < m — 1. So the error propagates more slowly 
than the exclusion of the states. 

The closed form solution for p m (i) reads 

i—m N ' N ' v 

and its derivation is presented in Appendix [A] An example is plotted in Fig. [3] (a). 

However, the expression above, although exact, is fairly complicated. Let's analyze the mean number of 
mafia players after t days and nights, that is 

M 

(m)(t) = J2mp m (t). (15) 

m— 

In birth and death processes with transition coefficients proportional to state's label (as radioactive decay 
process or Yule process) such mean may be obtained by a straightforward calculation with the use of (|13l) 

{ m)(t) = ———————— (to) (t - 1) = fn N ~ 2i ~ 1 ) M. (16) 

V /W N-2t + 2 K n ' 111 N-2i J y ' 

The result is correct only when N — 2t — M > 0, that is, the average is taken over probabilities of the correct 
states. 



<s 



B. Continuous time approximation 



Usually differential equations are simpler than difference equations. In this Section we crudely approximate 
the discrete-time evolution (fTo| by its continuous-time version. We know that 

n,f\ aa , p m (t + A) —pm(t) d 
Pm (t + 1)- Pm (t) 4^ -p m (t). (17) 

Let's change in (ff"3| the difference p m (t + I ) — p m (t) into the differential j^p m {t) 

s p»(*) = - jTTat^W + iv-^ Pm+l(i) (18) 

where m G N and with the initial condition p m (0) — 5 m M- When a relative change of the function derivative 
is over unit length interval [t, t + 1], the approximation should be good p m (t) p m (t). Intuitively speaking, 
when the chance of lynching a mafia member during one day is small, such an approximation is justified. 
However, in the course of this paper we do not estimate the approximation error. 



The solution of the differential equation (jT8j) reads 



. \ M—m / . \ m 



in 



see Appendix [B] for detailed derivation. This is more convenient to work with than the exact expression 
(fT4f . An example is plotted in Fig. [3] (a). 

For example, maxima of p m (t) are easily found 

N ( /m\2\ 



The mean number of mafia players, defined as in (|f 5[) . is 



(m)(t) = Myi-|, (21) 

which is a simpler expression than its analogue (|f 6p . The comparison of the discrete-time formula with its 
continuous-time approximation is in Fig. [3](b). 




- turn 



Figure 3: Dynamics of the Mafia game for the initial conditions N = 32 and M = 4. We compare results for the 
discrete-time evolution with its continuous-time counterparts, (a) The probability that after t turns there are exactly 
m mafia members alive — the exact results p m (t) (dots) and their approximations p m (t) (lines), (b) The mean 
number mafia members — the discrete result (m)(t) (dots) and its continuous time variant (rh)(t) (lines). 
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It is tempting to estimate the mafia winning-chance. As po (t) is the probability that the citizens won, setting 
time for which there is only one player left gives the total citizen winning-chance: 



w(n, m) = 1 — po (V) = 1 — ( 1 — \/l — ~~ — ~ I ( 22 ) 



rn 



Not surprisingly, in the continuous-time approximation, there is no explicit dependence on parity (the discrete 
properties are lost). 



VI. GENERAL SOLUTION 



The result for p m (t) (I14p is not only interesting by itself — it also gives us the solution for the mafia winning- 
chance Note that po(t) is the probability that after t turns citizens win. When total number of players is 
odd, game lasts for t dd = (N — l)/2 days and nights, leaving 1 player alive. Thus po(t dd) is the probability 
the citizen win the game. For an even number of players, after t even = N/2 turns there are no players alive. 
But now po(t even ) is the probability that there are no mafia members, so in fact — no-one was killed during 
the last night and citizen win. So w(n,m) = 1 — po( JV ~ w 2 mod 2 )- Consequently, the closed form formula for 
the mafia winning- chance reads 

w(n, m) = 1 - Y ( m ) (-If „,, (n ~ ^" r-77 . (23) 

v ' ' J n\\((n mod 2) - i)U 



Let's show what we have proven for a game with only one mafia member — the asymptotic behavior @. 

To get the asymptotic formula for w(n, to) it suffices to notice that in the limit n — ¥ oo (taken over a selected 
parity) only the first two terms of the sum contribute. For i = the term is equal to 1. The terms with 
i > 1 do not matter as (n — i)\\/(n — 1)!! — > 0. Consequently, we may write 

(n - 1)" 

w(n,m) « 771- r~ — (24) 



ni 

lC \ (n mod 2)-l/2 m 



where we used the result from game with one mafia member, that is, approximating §3§ with (0). The 
approximate formula can be formally expressed as 

It is tempting to ask for the optimal number of mafia members m op t, for a given number of players. An 
interesting game is one in which both groups have the same chance of winning 



l opt ~ 2 (2 J ^ *- 25 ^ 



which are plotted in Fig. |4] 



VII. CONCLUSION 



In this paper we make, and solve, a mathematical model of the Mafia game. The starting point of our 
calculations is random lynching process ([1]). 

The main interest of this research is analysis of the mafia winning-chance w(n,m) ([2]). We find the closed- 
form solution (|23|) . as well as its asymptotic approximation (|24jl . We prove some discrete properties of 
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™opt 




Figure 4: Optimal number of mafia members for a given number of players. Points show numerical results, that is, 
number of mafia members m op t for which w(n,m op t) is the closest to 1/2. Lines show the approximation (|25[) . Dots 
and solid line are for even number of players, whereas squares and dashed line — odd. 

w(n,m), namely that adding two citizens decreases the mafia winning-chance, whereas adding a citizen 
(which is an odd player) — increases. 

The dynamics of the game is also of our interest. We calculate the probability that after t days and nights 
there is exactly a certain number of mafia members alive (|13[). We obtain closed- form expressions (|14|) and 
their convenient approximations (|18|) . 

In every realistic Mafia game, there is a psychological part. A dry mathematical model may be a backbone 
of a more complex one, but certainly, is not enough to describe a real game. During course of play, citizens 
gain some information — either by discovering another's identity by themselves, or trying to catch messages. 
Furthermore, voting may be subject to some kind of witch-hunt mentality. Moreover, rarely are all players 
the same — usually there are ones with higher and lower influence on the others. 

All these processes may be investigated mathematically, or numerically. There is still a lot of challenge in 
finding an appropriate model and obtaining the results. Comparison of the theory with the experiment may 
be crucial. 

Out of author's personal experience, in games when players don't know each other, usually a random person 
is killed (e.g. a person who is always nervous, not just because of being in the mafia). When they know each 
other very well, very dim signals can be used to reveal one's identity. 
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Appendix A: Solution of difference equation for p m (t) 



We want to solve the following set of difference equations (|13l) . 

. , (N - 2t) - m , . m + 1 , , 
p m (t + l) = N _ — p m [t) + N _ 2t Pm+i{t) ■ m^TL 

with the initial condition p m (0) = S m M- Let's introduce a generating function 

oo 

F(t,z)= ^ Pm (t)z m , (Al) 

with the respective initial condition F(0,z) = z M . After side-by-side multiplication of ([13]) by z m and 
summation we get 

/ x z dFCt.z) 1 dFCt.z) ,„ , 

N-2tdz v ' ; 



B(t) 

The linear operator B(i) has the eigenvectors of the form (z — l) fc with the corresponding eigenvalues 
(1 — k/(N — 2i)). As the eigenvectors do not depend on t, the computation of F(t, z) is straightforward: 

F(t, z) = B(t - l)B(t - 2) • • • B(0)z M (A3) 

-f ^<*- 1 )'( 1 - 3f ^)(i-s^)"-('-s 



i=0 
JVf 



fc=0 



M\ ^(AT-2*)!! (JV-i)!! 



A!! (JV-2*-i)H 



Applying once more the binomial theorem, and comparing the result (IA3[) to the definition of the generating 
function (jAll) we get 



' Nil (N-2t-i)\\' 



or flUl. 



Appendix B: Solution of differential equation for p m (t) 

We want to solve the following set of differential equations (fT5|) . 

d m m + 1 

with the initial condition p(0) = S m j\j. Once again, we employ a generating function 

oo 

G{t,z)=Y J Prn{t)z m , (Bl) 



m=0 

M 



with the respective initial condition G(0, z) = z . We get the partial differential equation 

dG(t,z) _ (-£_+!) dGj^z) 
di ~ N -2t dz ' 



(B2) 
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which we want to solve by the method of characteristics. Let's have integral curves of G(t, z) in the form of 



0G(t, z) dt dG(t, z) dz 
dt dip dz dip 



0. 



Comparison of (|B2[) with (|B3|) yields in 



(B3) 



dt dz 

dip 



N-2t 



z-V 



(B4) 



or 



G(t, z) = f (in VN - 2t(l - z) ) 



(B5) 



All we need is to find a function f(x) that satisfies the initial condition. The explicit form of the generating 
function reads 



G(t, Z )= (l-y/T= ?(l-2)) 



M 



(B6) 



which gives the result 



M-r 



_ 2t 



or 0. 
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